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We explore the non-linear dynamics of two parallel periodically patterned metal surfaces that 
are coupled by the zero-point fluctuations of the electromagnetic field between them. The resulting 
Casimir force generates for asymmetric patterns with a time-periodically driven surface-to-surface 
distance a ratchet effect, allowing for directed lateral motion of the surfaces in sizeable parameter 
ranges. It is crucial to take into account inertia effects and hence chaotic dynamics which are 
described by Langevin dynamics. Multiple velocity reversals occur as a function of driving, mean 
surface distance, and effective damping. These transport properties are shown to be stable against 
weak ambient noise. 
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The observation of Casimir forces in increasingly small 
devices on sub-micron scales has generated great cur- 
rent interest in exploring the role of these forces for 
the development and optimization of micro- and nano- 
electromechanical systems [H, 0, H|. These systems can 
serve as on-chip fully integrated sensors and actuators 
with a growing number of applications. It was pointed 
out that Casimir forces can make an important contri- 
bution to the principal cause of malfunctions of these 
devices in form of stiction that results in permanent ad- 
hesion of nearby surface elements Q ■ This initiated in- 
terest in repulsive Casimir forces by modifying material 
properties as well as the geometry of the interacting com- 
ponents d, @. 

A complementary strategy, which will be considered 
in this Letter, is to make actually good use of Casimir 
forces between metallic surfaces and employ them to ac- 
tuate components of small devices without contact. We 
will demonstrate that this can be achieved by coupling 
two periodically structured (locally) parallel surfaces by 
the zero-point fluctuations of the electromagnetic field 
between them. The broken translation symmetry paral- 
lel to the surfaces results in a sideways force which has 
been predicted theoretically 0, [1| and observed experi- 
mentally between static surfaces Q. If at least one of 
the surfaces is structured asymmetrically there is an ad- 
ditional breaking of reflection symmetry and the surfaces 
can in principle be set into relative lateral motion in the 
direction of broken symmetry. The energy for this trans- 
port has to be pumped into the system by external driv- 
ing. This can be realized by setting the surfaces into 
relative oscillatory motion so that their normal distance 
is an unbiased periodic function of time. Since the side- 
ways Casimir force decays exponentially with the normal 
distance Q, the surfaces experience an asymmetric peri- 
odic potential that varies strongly in time. 

This scenario resembles so-called ratchet systems 
that have been studied extensively du ring the last decade 
in the context of Brownian particles , molecular mo- 



tors [T3 | and vortex physics in superconductors (l3| . to 
name a few recent examples. Most of the works on ratch- 
ets consider an external time-dependent driving force act- 
ing on overdamped degrees of freedom to rectify thermal 
noise. For nano-systems, however, it has been pointed 
out that inertia terms due to finite mass should not be 
neglected and, actually, can help the ratchets to perform 
more efficiently than their overdamped companions [l4l ]. 
Finite inertia typically induce in Langevin dynamics de- 
terministic chaos that has been shown to be able to mimic 
the role of noise and hence to generate directed trans- 
port in the absence of external noise [la ]. Here we use 
this effect in the different context of so-called pulsating 
(or effectively on-off ) ratchets where the strengths of the 
periodic potential varies in time [Io| . We consider weak 
thermal noise only to test for stability of the inertia in- 
duced transport — not as the source of driving. We note 
that in the absence of inertia, finite thermal noise is nec- 
essary for on-off ratchets to generate directed motion. 

In this Letter we demonstrate that the system de- 
scribed above indeed allows for directed relative motion 
of the surfaces due to chaotic dynamics caused by the 
lateral Casimir force. The transport velocity is stable 
across sizeable intervals of the amplitude and frequency 
of surface distance oscillations and damping. The veloc- 
ity scales linear with frequency across these intervals and 
is almost constant below a critical mean distance beyond 
which it drops sharply. The system exhibits multiple 
current reversals as function of the oscillation amplitude, 
mean distance and damping. The Casimir ratchet allows 
contact-less transmission of motion which is important 
since traditional lubrication is not applicable in nano- 
devices. This actuation mechanism should be compared 
to other actuation schemes as magnetomotive or capaci- 
tive (electrostatic) force transmission. The Casimir effect 
induced actuation has the advantage of working also for 
insulators and does not require any electrical contacts 
and/or external fields. It can also scale down effects of 
parasitic capacitance that reduces the efficiency of actua- 
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tion at high frequencies [4|. It should be mentioned that 
other applications of zero-point fluctuation induced (van 
der Waals) interactions to nano-devices have been exper- 
imentally realized already to construct ultra-low friction 
bearings from multiwall carbon nanotubes (ljj . 

In the following, we consider two (on average) paral- 
lel metallic surfaces with periodic, uni-axial corrugations 
(along the xi-axis) that have distance H, see inset (a) of 
Fig. Q] To begin with, we assume that both surfaces are 
at rest with a relative lateral displacement b. Then the 
surface profiles can be parametrized as 

oc 

hiipn) = aJ2 c ne 27rmxl/Xl +c.c, (la) 

n=l 
oc 

h 2 (x 1 ) = aJ2 d ne 27rm(xi ~ b)/X2 +c.c, (lb) 

n=l 

where a is the corrugation amplitude, Ai, A2 are the cor- 
rugation wave lengths, c n , d n arc Fourier coefficients. 

The Casimir energy £ of this configuration is the 
change of the ground state energy of the electromagnetic 
field due to the suppression of the tangential electric field 
at the surfaces. The dependence of £ on H and b causes 
macroscopic forces on the surfaces. For a varying sepa- 
ration H this is the normal Casimir attraction between 
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metallic surfaces |17| . modified by the corrugations. Be- 
low we will assume H = H(t) to be a time-dependent dis- 
tance that is kept at a fixed oscillation by an additional 
external force from clamping to an oscillator. In such 
setup, the surfaces can react freely only to the lateral 
force component T\ a t(b, H) = —d£/db. This side- ways 
force has been computed for sinusoidal corrugations to 
second order in the amplitude a, using a path integral 
formulation This result is readily extended to peri- 
odic profiles of arbitrary shapes as described by Eq. J!} . 
We find that the corrugation length have to be commen- 
surate, A1/A2 = p/q with integers p, q in order to produce 
a finite lateral force per surface area. For the purpose of 
this work, it is sufficient to consider the case p = 1 . The 
lateral (6-dcpcndcnt) part of the Casimir energy per sur- 
face area can then be written as 



£{b) 
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to order a 2 . The exact form of the function J was ob- 
tained in Ref. [H in terms of transcendental functions. For 
the present purpose it is sufficient to use the simplified 
expression 



J(u) 
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(3) 



with 7 = 12.4133, which is exact for both asymptotically 
large and small u and approximates the exact results 
with sufficient accuracy for all u. (The maximal devi- 
ation from the exact result is « ±0.5% around u = 0.5.) 
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FIG. 1: The lateral Casimir force acting between the two 
surfaces as function of the shift b at time s = and half 
period s = tt/uj (drawn to a larger scale by a factor 10 3 ) 
for parameters n = 0.65, Ho = 0.1A. Insets: (a) Surface 
profiles at their equilibrium position at b = 0.182 (b) Periodic 
variation of the maximum force at b — with time. 

The Casimir potential of Eq. @ has two interesting prop- 
erties which are useful to the construction of a ratchet. 
First, it decays exponentially with H, and thus can be 
essentially switched on and off periodically in time by 
oscillating H. Second, the potential is not only periodic 
in b but acquires asymmetry from the surface profiles 
at small H <C A and an universal symmetric shape for 
H ^> A since the effect of higher harmonics of the surface 
profile is exponentially diminished (l8j . 

The relative surface displacement b(t) can be consid- 
ered as a classical degree of freedom with inertia. Its 
equation of motion is described by Langevin dynamics of 
the form 



pb + 1P b = Tut[b, H(t)] + V^prat) ■ 



(4) 



where p is the mass per surface area, 7 the friction co- 
efficient, T the intensity (divided by surface area) of 
the Gaussian noise £(i) with zero mean and correla- 
tions (£(t)£(^')) = S(t — t') so that the Einstein rela- 
tion is obeyed. This stochastic term describes ambient 
noise due to effects of temperature and pressure. (Ad- 
ditional contributions from thermally excited photons to 
the Casimir force can be neglected at surface distances 
well below the thermal wavelength Hc/(2T).) The sys- 
tem is driven by rigid oscillations of one surface so that 
the distance H(t) = Hog(t) oscillates about the mean 
distance Hq with g(t) = 1 — ?]cos(fW). For simplicity, 
we consider now equal corrugation lengths Ai = A2 = A. 
We define the following dimensionless variables: b = b/X, 
s = t/r for lateral lengths and time with the typical time 
scale t = {\/a)\J ' pH^/hc resulting from a balance be- 
tween inertia and Casimir force. Hence velocities will 
be measured in units of vq = A/r. There are five di- 
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mcnsionless parameters which can be varied indepen- 
dently for fixed surface profiles: the damping 7 = 7*7, 
the angular frequency to = rO, the driving amplitude 
77, the scaled mean distance H /X and the noise inten- 
sity T = (T ' I 'Hc)(Hq I a 2 ) . The dimensionless equation of 
motion for b(s) reads 



b + jb = T lat [b,g(s)} + J2jTas) 



with the Casimir force 

47T 



^Utt(b, 9) = — ^2 f™ C0S ( 2lT nb)J f ng 



fJ 
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(5) 



(6) 



where we have chosen surface profiles with c n = 
iy/ f n / (2fi), d„ = \J f n / (2n) with real coefficients /„ in 
Eq. fl}, and #(s) = 1 — 77COs(ws). 

We begin to analyze Eq. ([5]) by noting that directed 
transport is possible in certain parameter ranges even in 
the deterministic case where noise is absent. However, 
to probe the robustness of transport, we consider in the 
following primarily the limit of weak noise by choosing 
T = 1CP 3 . In fact, it has been shown for undcrdamped 
ratchets with time-independent potentials and periodic 
driving that even an infinitesimal amount of noise can 
change the rectification from chaotic to stable 



look for similar generic behavior of our pulsating ratchet, 
we consider a specific geometry consisting of a symmetric 
and a sawtooth-like surface profile corresponding to three 
harmonics with f x = 0.0492, f 2 = 0.0241, f 3 = 0.0059 
and f n = for n > 3. Inset (a) of Fig. Q] shows these 
profiles in their stable position with b = 0.182 that mini- 
mizes the Casimir energy The resulting spatial variation 
of the Casimir force with b is plotted in Fig.[T]for minimal 
(s = 0) and maximal (s = tt/lu) surface distance with pa- 
rameters Ho/A = 0.1, 77 = 0.65. It can be clearly seen 
that the asymmetry is reduced at larger distance where 
the variation of the force becomes more sinusoidal. In- 
set (b) shows the on-off-likc time-dependence of the force 
amplitude at b = due to the oscillating surface distance. 

The non-linear equation of motion of Eq. |j5J) has to 
be solved numerically. The trajectory b(s) was obtained 
from a second order Rungc-Kutta algorithm. As initial 
conditions we used an equidistant distribution over the 

interval [—1,1] for b(0) and 6(0) = 0. For each set of 
parameters we calculated 200 different trajectories from 
varying initial conditions and noise, each evolving over 
4 x 10 3 periods 2ir/oj so that transients have decayed. 
The average velocity ((v)} involves two different averages 

of b(s): The first average is over initial conditions and 
noise for every time step, then the averaged trajectory is 
averaged over all discrete times of the numerical solution. 
For an efficient directed transport it is not sufficient to 
have only a finite average ((f)). To exclude trajectories 
with a high number of velocity reversals, the fluctuations 




FIG. 2: Mean {(v)) and standard deviation a of the (nega- 
tive) velocity as function of the driving amplitude r\ for the 
frequencies u) — 5.0 and uj — 4.72 (for the latter only the 
stable plateau is shown). The parameters are Ho = 0.1 A, 
7 = 0.9, T — 10~ 3 . Inset: Dependence of the same quan- 
tities on frequency for fixed 77 = 0.65. Straight dashed lines 
correspond in both graphs to the velocity ui/(2n). 



1M- To about the average velocity must be small, i.e., the vari- 
ance cr 2 = ((v 2 )} — ((v)} 2 must be smaller than ((v)) 2 . 

Naively, one can expect directed motion of the sur- 
face profile h,2(xi) into the positive 2; 1 -direction (6 < 0) 
since the Casimir force in Fig. [T] is asymmetric with neg- 
ative values lasting for longer time than positive ones. 
However, the actual behavior is more complicated due to 
chaotic dynamics. Fig.[2]shows the dependence of the av- 
erage velocity and its standard deviation a on the driving 
amplitude 77 and frequency to for Ho = 0.1 A, 7 = 0.9. For 
a fixed frequency there is an optimal interval of driving 
amplitudes across which the average velocity is almost 
constant with ((«}} ~ —u)/(2tt). Small deviations from 
the latter value result from noise as we have checked by 
studying the dynamics at T = 0. At higher driving am- 
plitudes we observe a second narrower interval with max- 
imal ((d)) which is more strongly reduced and smeared 
out from its deterministic value —2 x lu/(2tt) by noise. 
At the plateaus of constant velocity the standard devia- 
tion a is substantially reduced, rendering transport effi- 
cient. Outside the plateaus velocity reversals occur and 
tr increases linearly with 77. For fixed amplitude 77, the 
average velocity is stable at the value — lj/(2it) over a 
sizeable frequency range (see inset of Fig. [J) . 

In order to understand the observed behavior we have 
analyzed the dynamics in the three dimensional extended 
phase space. There attractors of the long-time dynamics 
can be identified from Poincare sections using the period 
2tt/lo of the surface oscillation as stroboscopic time. To 
obtain a compact section, the trajectory is folded period- 
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FIG. 3: Mean {{v)) and standard deviation <r of the (negative) 
velocity as function of (a) the mean plate distance Ho for 7 = 
0.9 and (b) damping 7 for Ho — 0.1A. The other parameters 
are 77 = 0.65, u = 5.0, f = 10 -3 . 

ically in x\ on one period of the Casimir potential. From 
these sections we can distinguish between periodic and 
chaotic orbits. As a start, we consider the deterministic 
limit with T — 0. The plateaus around 77 = 0.65 and 
i] = 0.7 result both from periodic orbits of period one, 
corresponding to a single point in the Poincare section. 
On the right (downward) edges of the first plateaus we 
have observed period doubling, i.e., a periodic attractor 
with period two. Upon a further increase of rj, chaotic 
orbits dominate the motion. Hence the system exhibits 
a period-doubling route to chaos with enhanced velocity 
fluctuations. The findings apply basically also to weak 
noise (T = 10~ 3 ) but the sharp points of the periodic at- 
tractors in the Poincare sections are smeared out leading 
to a decreased ((v)) . The transition from chaotic to peri- 
odic dynamics at the beginning of the rising edge of the 
plateaus is accompanied by a velocity reversal. This is 
consistent with the earlier observation for non-pulsating 
potentials that velocity reversals are due to a bifurcation 
from chaotic to periodic dynamics (l5| . 

The amplitude of the Casimir potential can be tuned 
by varying the mean distance Ho. From Fig.[3la) we see 
that the dynamics show a sharp transition at a critical 
Hq/X from efficient transport with large ((v)) and small 
a to chaotic dynamics with vanishing velocity. The tran- 
sition is accompanied by a velocity reversal and peaked 
velocity fluctuations. Interestingly, below the transition 
((v)} is almost constant independently of Hq/X. The ob- 
served transport behavior is also stable against a change 
of effective damping 7 as shown in Fig. [3](b) . Whereas 
fluctuations increase with decreasing 7, there is a stable 
plateau of constant average velocity across which fluc- 
tuations are diminished. In the deterministic limit, we 
have also observed additional plateaus with inverted and 
doubled average velocity by varying 7 and r\. Remnants 
of a second plateau around 7 = 1.9, washed out by noise, 
can be seen in Fig. [3jb). 

Finally, let us estimate typical velocities vq = X/t. 
With the typical lengths H$ = 0.1/im, a = lOnm real- 
ized in recent Casimir force measurements [9J and an 



area mass density of p = 10g/m 2 for silicon plates 
with a thickness of a few microns, one obtains vq = 
-J hca 2 1 P-ffp ~ 5.5mm/s. The actual average velocity 
Vqw/2-k is of the same order for the frequencies studied 
above. For A = l/_tm, the time scale is r = A/u ~ 10~ 4 s 
leading to driving frequencies and damping rates in the 
kHz range for the parameters considered here. 

Our results show that Casimir interactions offer novel 
contact-less translational actuation schemes for nanomc- 
chanical systems. Similar ratchet-like effects are ex- 
pected between objects of different shapes as, e.g., peri- 
odically structured cylinders, inducing rotational motion. 
The use of fluctuation forces appear also promising to 
move nano-sized objects immersed in a liquid where elec- 
trostatic actuation is not possible. Another application 
is the separation and detection of particles of differing 
mass adsorbed to the surfaces. For surfaces oscillating at 
very high frequencies additional interesting phenomena 
related to the dynamical Casimir effect occur [l9|], lead- 
ing to the emission of photons that could contribute to 
ratchet-like effects as well. 

I acknowledge useful discussions with R. Golestanian 
and M. Kardar. 
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